Solar neutrinos and the solar composition problem 
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Standard solar models (SSM) are facing nowadays a new puzzle: the solar composition problem. 
New determinations of solar metal abundances lead SSM calculations to conflict with helioseismo- 
logical measurements, showing discrepancies that extend from the convection zone to the solar core 
and can not be easily assigned to deficiencies in the modelling of the solar convection zone. We 
present updated solar neutrino fluxes and uncertainties for two SSM with high (old) and low (new) 
solar metallicity determinations. The uncertainties in iron and carbon abundances are the largest 
contribution to the uncertainties of the solar neutrino fluxes. The uncertainty on the ^''N(p,7)^^0 
rate is the largest of the non-composition uncertainties to the CNO neutrino fluxes. We propose 
an independent method to help identify which SSM is the correct one. Present neutrino data can 
not distinguish the solar neutrino predictions of both models but ongoing measurements can help 
to solve the puzzle. 



This paper is part of a series led by John Bahcall that 
spans more than 40 years The goals of this series 
have been to provide increasingly more precise theoreti- 
cal calculations of the solar neutrino fluxes and detection 
rates and to make increasingly more comprehensive eval- 
uations of the uncertainties in the predictions. 

In the past, standard solar models (SSM) had been 
steadily showing a discrepancy with neutrino measure- 
ments, named the solar neutrino problem, that defini- 
tively led to the discovery of flavor change in the lepton 
sector. The significant difference of the charged current, 
electron scattering and neutral current measurements of 
the epochal Sudbury Neutrino Observatory (SNO) and 
Super-Kamiokande experiments have demonstrated the 
new phenomenon. As soon as the solar neutrino puzzle 
was solved, a new and different puzzle has emerged. 

Recent refined determinations of the abundances of 
heavy elements on the solar surface have led to lower 
heavy-element abundances 0] than previously measured. 
Solar models that incorporate these lower abundances 
conflict dramatically with helioseismological measure- 
ments. Detailed studies of the physical inputs have been 
done to resolve this controversy; with particular focus on 
the region where more dramatically low metallicity solar 
models fail to explain helioseismological measurements, 
right below the solar convective envelope. But, recent 
analysis of precise low-degree acoustic oscillations of the 
Sun measured by the Birmingham Solar-Oscillations Net- 
work (BiSON) indicate that the discrepancies between 
solar models constructed with low metallicity and helio- 
seismic observations extend to the solar core and thus 
cannot be attributed to deficiencies in the modelling of 
the solar convection zone [3]. 

A number of studies have been done to refine the esti- 
mated uncertainties and search for systematic effects that 
could account for the disagreement of helioseismological 
predictions of the low metallicity solar models. Today, 



we lack a solution to this puzzle and extra handles are 
needed to identify the source of the discrepancy. In this 
letter, we propose a new method to discriminate whether 
the low heavy-element abundances should be accepted 
with the other solar inputs in the SSM or not. 

Here we present the most updated solar model calcula- 
tions with the two solar abundance determinations, high 
metalhcity labelled as BPS08(GS) [1] and low metallicity 
labelled as BPS08(AGS) ^. There are two steps forward 
that contribute to major improvements in the estimated 
precision of the neutrino flirxes: i) improved accuracy 
of the •^He-^He cross section and ii) reduced systematic 
uncertainties in the determination of the surface compo- 
sition of the Sun. 

Previously, the uncertainty on the S-factor of the 
■^IIe(Q;,7)^Be reaction was due to an average discrepancy 
in results obtained by the detection of the delayed 7 rays 
from ^Be decay and the measurement of prompt 7 emis- 
sion. LUNA collaboration reported on a new high pre- 
cision experiment using both techniques simultaneously 
at the lowest interaction energy ever reached in the lab. 
The S-factors from the two methods do not show any 
discrepancy within the experimental errors, obtaining 
S(0)=0.567±0.018±0.004 keV barn [|. With this new 
determination, the uncertainty on the predicted ^B (^Be) 
neutrino flux due to S'34 is reduced from 7.5% to 2.7% 
(8.0% to 2.8%). Best-estimates and l-cr uncertainties of 
the other important astrophysical factors are summarized 
in Table 1 of Ref . @ , with minor modifications J] . 

Previously, we preferred solar model calculations with 
conservative uncertainties in the element abundances, es- 
timated by the difference between the high and low metal 
abundances, although it was recognized that estimated 
uncertainties based on each analysis independently were 
smaller [1]. New detailed studies of helioseismological 
data [3] have shown that low-degree helioseismology data 
clearly allow to discriminate between predictions for the 
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TABLE I: Predicted solar neutrino fiuxes from solar mod- 
els. The table presents the predicted fluxes, in units of 
10^°(pp), lO^(^Be), 10**(pep,i3N,i5o), 10''(*B/^F), and 
10^ (hep) cm-^s"\ Columns 2 and 3 show BPS08 for high 
and low metalicities; and column 4 the flux differences be- 
tween the models. 



Source BPS08(GS) BPS08(AGS) Difference 



pp 


5.97(1 ± 0.006) 


6.04(1 ± 0.005) 


1.2% 


pep 


1.41(1 ± 0.011) 


1.45(1 ± 0.010) 


2.8% 


hep 


7.90(1 ± 0.15) 


8.22(1 ± 0.15) 


4.1% 


^Be 


5.07(1 ± 0.06) 


4.55(1 ± 0.06) 


10% 


«B 


5.94((1 ± 0.11) 


4.72(1 ± 0.11) 


21% 


13N 


2.88(1 ± 0.15) 


1.89(1 tVit) 


34% 


15q 


2.15(1 


1.34(1 t°-\l) 


31% 


"F 


5.82(1 


3.25(1 


44% 



CI 


8.4610 


B7 

88 


O.»D_0.70 


Ga 


127.91^ 


.1 
.2 


120. 5i?-^ 



solar core structure from solar models with high and low 
metallicity. High metallicity models agree with helioseis- 
mological data but not with recent determinations of so- 
lar abundances. On the contrary, low metallicity models 
are inconsistent with helioseismological data but use the 
recent determinations of solar abundances. That is, high 
and low metallicity solar models really are two different 
groups of models, and conservative uncertainties defined 
as described above are, in fact, too conservative. We find 
that uncertainties have to be estimated for each class in- 
dependently. In what follows, we adopt for each solar 
abundances set 0, 13| both the recommended values and 
uncertainties for each chemical element. 

Table H] presents, in the second and third columns, 
labelled BPS08(GS) and BPS08(AGS), our best solar 
model calculations for the neutrino fluxes. Both mod- 
els were calculated with the new input data for nu- 
clear physics 0, Q. BPS08(GS) was computed with 
the high (old) metallicity solar abundances determina- 
tions '4| and BPS08(AGS) with the most recent analyses 
of solar abundances 2-| that lead to conflicts with helio- 
seismological measurements. 

Table U also shows the overall uncertainties in the neu- 
trino fluxes. Uncertainties from different sources have 
been calculated using power-law expansions with coeffi- 
cients given in Table |TT] and added quadratically. Un- 
certainties due to errors in the solar composition (met- 
als) have been determined for each relevant metal sep- 
arately 01 and abundance uncertainties are taken from 
the respective solar composition compilations @] . This is 
a more consistent way for estimating fiuxes uncertainties 
and leads to, incidentally, reduced overall uncertainties. 
This result is particularly relevant in the context of CNO 
fluxes and the solar abundance problem. Solar composi- 
tion dominates the uncertainty in the CNO fluxes (13% 
for i^N, 12% for and 17% for "F) and it is not dom- 



inant for all other fluxes. Among the non-composition 
sources, the most important contribution to the total er- 
ror in the CNO fluxes is due to the 5*1, 14 uncertainty. 

contributes with 6% (8%), while diffusion and opac- 
ity contribute with 4% and 5% (5% and 6%) to the ^^N 
(^^O) neutrino flux. For the ''Be flux, opacity (3.2%) con- 
tributes most to the total error, followed by S34 and 6*33 
(2.8% and 2.5%); while for ^B, opacity (6.8%), diffusion 
(4.2%) and Sn (3.8%) dominate the uncertainties. 

A careful look to the logarithmic partial derivatives 
with respect to the individual solar abundances and the 
uncertainties of the individual elements shows that the 
major contribution to the uncertainties in the ^"^N and 
^^O neutrino fluxes largely comes from the uncertainty 
in the carbon abundance. Therefore, a reduction of the 
carbon abundance uncertainty directly translates into a 
substantial improvement of the CNO neutrino flux uncer- 
tainties. On the other hand, ^B and ^Be fluxes are much 
more sensitive to the iron abundance (due to its influ- 
ence on the core opacity) and therefore a reduction of 
the iron abundance uncertainty directly translates into a 
substantial improvement of the *B and ^Be neutrino flux 
uncertainties. 

Finally, in Table |TT] we also included power-law coeffi- 
cients for the dependence on two important helioseismic 
quantities on solar model input parameters, which show 
the relevant inputs that contribute to their uncertain- 
ties. SSM predictions for these quantities are Rcz — 
0.713(0.728) ± 0.004 and Ys = 0.243(0.229) ± 0.035 for 
the high (low) metallicity models; while helioseismolog- 
ically determined values are Rcz — 0.713 ± 0.001 and 
Ys = 0.2485 ±0.0034 0. 

Opacity uncertainties have been conservatively com- 
puted by defining the 1-a uncertainty as given by the 
flux differences between two SSMs that have been com- 
puted one with standard opacities and the other with 
opacities in the solar interior increased by 2.5%. Such a 
difference in the solar core opacities is representative of 
the maximum difference between the two most up-to-date 
radiative opacitity calculations, OPAL and OP 



Hi 



Based on the results presented above, we make two 
recommendations for future work to improve the uncer- 
tainties in the SSM neutrino fluxes. First, the uncer- 
tainty in the low energy extrapolation of the rate of the 
^^N(p,7)^^0 reaction should be reduced to below 5% (al- 
most a factor 2 improvement) in order that the uncer- 
tainty due to astrophysical cross sections does not dom- 
inate the non composition uncertainties in any of the 
neutrino fluxes. Second, the uncertainty in the surface 
element abundances iron and carbon should be reduced 
to ± 0.02 dex (factor 2 to 3 improvement) to significantly 
improve the radiative opacity and the CNO composition 
uncertainties. The requirements are tough, but the out- 
come is rewarding. 

Table |T] shows that the central values of the neutrino 
fluxes in the two models BPS08(GS) and BPS08(AGS) 
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TABLE II: Logarithmic partial derivatives a{i,j) of neutrino fluxes and helioseismic quantities with respect to nuclear cross 
sections, solar luminosity and age, element diffusion rate, radiative opacity and solar abundances. 



Source Sii S33 S34 Si7 Shop Si, 14 St^^,, Lq Age Diff Opac C N O Ne Mg Si S Ar Fe 

pp 0.090 0.029 -0.059 0.000 0.000 -0.004 0.000 0.808 -0.067 -0.011 -0.099 -0.005 -0.001 -0.005 -0.004 -0.004 -0.009 -0.006 -0.001 -0.016 

pep -0.236 0.043 -0.086 0.000 0.000 -0.007 0.000 1.041 0.017 -0.016 -0.300 -0.009 -0.002 -0.006 -0.003 -0.002 -0.012 -0.014 -0.003 -0.054 

hep -0.112 -0.459 -0.072 0.000 1.000 -0.004 0.000 0.174 -0.118 -0.037 -0.398 -0.007 -0.002 -0.020 -0.014 -0.017 -0.036 -0.028 -0.005 -0.064 

''Be -1.07 -0.441 0.878 0.000 0.000 -0.001 1.000 3.558 0.786 0.136 1.267 0.004 0.002 0.053 0.044 0.057 0.116 0.083 0.014 0.217 

**B -2.73 -0.427 0.846 1.000 0.000 0.005 0.000 7.130 1.380 0.280 2.702 0.025 0.007 0.111 0.083 0.106 0.211 0.151 0.027 0.510 

i^AT -2.09 0.025 -0.053 0.000 0.000 0.711 0.000 4.400 0.855 0.340 1.433 0.861 0.148 0.047 0.035 0.051 0.109 0.083 0.015 0.262 

^^O -2.95 0.018 -0.041 0.000 0.000 1.000 0.000 6.005 1.338 0.394 2.060 0.810 0.207 0.075 0.055 0.076 0.158 0.117 0.021 0.386 

^''F -3.14 0.015 -0.037 0.000 0.000 0.005 0.000 6.510 1.451 0.417 2.270 0.024 0.005 1.083 0.061 0.084 0.174 0.128 0.023 0.428 



Rcz 


-0.061 


0.002 


-0.003 


0.000 


0.000 


0.000 


0.000 


-0.016 


-0.081 


-0.018 


-0.012 


-0.006 


-0.005 


-0.028 


-0.012 


-0.005 


0.002 


0.004 


0.001 


-0.009 


Ys 


0.134 


-0.005 


0.009 


0.000 


0.000 


0.001 


0.000 


0.373 


-0.110 


-0.073 


0.646 


-0.009 


-0.001 


0.023 


0.033 


0.037 


0.070 


0.048 


0.009 


0.089 



differ by 2-3 times the theoretical error expected within 
one model. This feature must be explored to estimate by 
how much solar neutrino experiments can contribute to 
solve the solar composition puzzle. Two of the neutrino 
fluxes (^B and ^Be neutrinos) have been measured with 
good precision, and this can lead to preliminary tests of 
the solar composition. SNO collaboration [l^i M\ 
has determined the *B neutrino flux with good precision 
by the measurement of the charged current and the neu- 
tral current detection of solar neutrinos in three phases of 
the experiment: I) heavy water target, II) heavy water + 
salt target and III) heavy water + ^He targets. All three 
phases lead to a very precise measurement of the ^B flux 
6%). Moreover, Super-Kamiokande electron scatter- 
ing data, when combined with the other solar neutrino 
data and the assumption of neutrino oscillations, con- 
tributes to further reduce the uncertainty in the *B flux. 
Since May 2007, the Borexino detector is collecting data 
on the electron scattering of ^Be neutrinos. The collabo- 
ration presented the latest results on ^Be solar neutrinos 
(^10% precision) based on 192 live days of data taking 
[l5|. With more statistics and better calibrations, the 
collaboration anticipates a final precision at the 5% level. 

In order to quantitatively determine how well the two 
solar models fit the neutrino data we proceed in two 
steps: a) use all of the available solar and reactor neu- 
trino data to determine the current constraints on neu- 
trino oscillation parameters and solar neutrino fluxes 
with the method described in Ref. b) use the so- 
lar neutrino fluxes inferred from data to test SSM fluxes. 
We discuss now which data have been included in the 
analysis and refer to Ref. for further details. So- 
lar neutrino data includes the average rate measured 
by the radiochemical experiments, chlorine 17| and gal- 
lium [3], the SuperKamiokande I zenith spectrum [3], 
the SuperKamiokande II day- night spectrum the 
SNO pure D2O phase day-night spectrum 12l |. the SNO 
salt phase day- night spectrum [l3|, the electron scat- 



tering, charged current and neutral current rates mea- 
sured in phase 3 IJ] and the Borexino ^Be and ^B 
rates [T^. Non solar neutrino data include the an- 
tineutrino data from KamLAND ^] and the marginal- 
ized x^(6'i3) function derived from the analysis of at- 
mospheric-|-K2K-f MINOS-fCHOOZ data The free 

parameters in the analysis are the neutrino parameters 
Am22, ^12 and ^13 and the reduced neutrino fluxes /b, 
/bc, /pp, /n and /o (/i is the i-th neutrino flux nor- 
malized to the BPS08(GS) prediction) subject to the lu- 
minosity constraint, i.e., the total luminosity produced 
by nuclear reactions (the p-p chains and the CNO cycle) 
equals the observed solar luminosity. We find that: 



tan^ 6*12 
/b 
/pp 



(7.7 ±0.2) X 10- 



= 0.46 



H-0.04 



0.05 

0.91 ±0.03 



= 1.00: 



-0.01 



,2 



sm t7i3 
/be 

LcNO 



6I13 = 0.014 



+0.011 

-0.009 

1.02 ±0.10 

+2.90 



= o.ol-% 



where Lqno is the fraction of solar luminosity due to the 
CNO nuclear fusion reactions. The oscillation parame- 
ters, central values and errors, are in excellent agreement 
with neutrino analysis where the solar neutrino fluxes 
are not treated as free variables 22, [2^. Compared to 
previous analysis [lit , the major improvement is in the 
determination of the ^Be flux. The Borexino measure- 
ments lead to a factor five improvement of the uncer- 
tainty in this flux. The upper bound 011 the CNO lu- 
minosity is not significantly different 



24l | because the 



improvement due to the ''Be determination leading to 
the more constrained ^^N and ^^O fluxes, compensates 
the lower CNO fluxes predicted by the solar model cal- 
culations using the most accurate cross section of the 
"N(p,7)i50 reaction used in BPS08. The results 
of this analysis allow us to experimentally determine 
the ratio of the termination chains in the p-p burning: 
R=(3He+4He)/(3He+3He)=0.19±0.02. 

The global analysis of neutrino data used above was 
discussed in Refs. [l^, HJ] with the goal to derive neu- 
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trino oscillation parameters and solar neutrinos fluxes 
from neutrino data independently of the solar model. We 
can now move to the second step in the analysis and 
test SSM predictions without dealing with solar neutrino 
data, which is encoded in the results presented above. In 
order to statistically test the hypothesis that there is no 
difference between a given solar model and the data, we 
define the function 

ij 

2 cx'p \ih 

<^ij ~ <^cxp,i<^cxp,jPij + <^th,iO'thjPij 

where aexp,i and ath,i are the experimental and theo- 
retical errors of the neutrino flux fi and p^J^ and pjj 
are the experimental and theoretical correlations of the 
fluxes i and j. Given the large uncertainties in some of 
the experimentally determined fluxes, we only include in 
the analysis the most precise ones (^B and ^Be fluxes), 
which have negligible experimental correlations. Theo- 
retical correlation coefficients for neutrino fluxes can be 
found in Ref. 8]. The analysis leads to: 

xIpso8{gs) = 0-9 (63%) and xIpsos{ags) = 1-5(47%), 

where the goodness of fit is shown in parenthesis and 
is used as a diagnostic for the fit of the model to the 
data. Neutrino data are not precise enough to reject the 
null hypothesis for any of the models, BPS08(GS) and 
BPS08(AGS). Adding the pp flux in the analysis does not 
improve the statistical power, changes to 0.95 and 1.6 
respectively. These results show that the theoretical and 
statistical errors are still too large to reject any of the two 
models. The analysis proposed here will guide us on how 
improvements can contribute to solve the composition 
problem. Eventually, the rejection of the null hypothesis 
would prompt us to use a likelihood ratio testing for a 
given parameter(s) connecting the models. More Super- 
Kamiokande data and the full analysis of the SNO data, 
with lower energy threshold, will sharpen the accuracy 
of the *B neutrino flux. Borexino ongoing measurements 
anticipate further improvements of the statistical and 
systematic errors which may help to discriminate mod- 
els. To illustrate what might be achieved, if we assume 
that future neutrino data will lead to /b ~ 1.0 ±0.02 and 
/bc = 1-00 ± 0.05, and theoretical errors on the *B and 
^Be fluxes are reduced from ±0.11 and ±0.06 to ±0.08 
and ±0.04 respectively, we obtain Xbpso8{gs) ~ ^ (^^ 
construction) and Xbpso8{ags) ~ ^-2 (4.5%), what will 
hint towards the relevance of doing an analysis of solar 
models as a function of the iron composition to all neu- 
trino data. Moreover, running and planned low energy 
experiments have a realistic chance to measure CNO neu- 
trinos and further distinguish the generic predictions of 
both models, given that predictions of BPS08(GS) and 
BPS08(AGS) for these fluxes differ about 30%. This kind 



of study, however, is beyond the scope of this paper. The 
hypothesis testing analysis with future data will show 
whether this effort is valuable to discriminate between 
the high and low composition SSM. 

In summary, standard solar models are facing the so- 
lar composition problem, i.e., new composition determi- 
nations, when included in solar models, conflict with he- 
lioseismological measurements. We are agnostic whether 
the new composition, the solar modelling, other inputs or 
a combination of them is the source of the conflict. We 
have proposed an independent method to identify which 
solar model is the correct one by using solar neutrino 
data. A method of analysis was proposed and prelimi- 
nary results have been derived. Ongoing solar neutrino 
measurements and forecast neutrino flux measurements 
will help in the search of the solution to this puzzle. 
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